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Note 


This report is a modified version of the paper (same title) 
presented at the Twelfth European Rotorcraft Forum, September 22-25, 
1986, Garmisch-Partenkirchen, Federal Republic of Germany. Errors in 
the blade section radius of gyration distribution used by the analysis 
have been corrected. Also, corrections to the measured pressure data 
have been made. 
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Abstract 


A comprehensive rotorcraft analysis model was used to predict 
blade aerodynamic and structural loads for comparison with flight test 
data. The data were obtained from an SA349/2 helicopter with an 
advanced geometry rotor. Sensitivity of the correlation to wake geom- 
etry, blade dynamics, and blade aerodynamic effects was investigated. 
Blade chordwise pressure coefficients were predicted for the blade 
transonic regimes using the model coupled with two finite-difference 
codes . 

SYMBOLS 


c.g. center of gravity 

C a blade section lift coefficient 

C n blade section normal force coefficient 

Cp blade chordwise surface pressure coefficient 

C'f/o ratio of rotor thrust coefficient to solidity 

r blade radial station 

R rotor radius 

r c /c ratio of tip vortex core size to blade chord 

a blade angle of attack 

a time derivative of blade angle of attack 

y advance ratio 

4> rotor blade azimuth station 


1. INTRODUCTION 


A joint rotorcraft research program has been established between 
the National Aeronautics and Space Administration (NASA) and the French 
Ministry of Defense. The purpose of the program is to correlate the 
NASA prediction code Comprehensive Analytical Model of Rotorcraft Aero- 
dynamics and Dynamics (CAMRAD) with flight test data from the French 
SA349/2 Helicopter and an advanced geometry rotor system. 

Structural and aerodynamic loads and performance were measured 
for several thrust levels and advance ratios. From the flight envelope, 
16 flight conditions were selected by NASA and Aerospatiale for documen- 
tation. In addition, data from three of these flight conditions were 
chosen for detailed correlation with CAMRAD. These conditions provided 
an opportunity to exercise various modeling features within CAMRAD. 

The pressure data from this flight test are a valuable addition 
to the limited data base of helicopter rotor blade pressures in forward 
flight, especially in the transonic regime. Predictions of shock forma- 
tions near blade tips using CAMRAD coupled with a finite-difference code 
have been restricted to comparisons with several model rotor data sets 
using rigid blades (described in Ref. 1). Chordwise pressure compari- 
sons presented here aid in further validating the coupling procedure. 

The flight test, the comprehensive model, the assumptions and 
approximations in modeling the aircraft and rotor in CAMRAD, and the 
resulting theoretical predictions compared with the data are described 
in this paper. Finally, recommendations for improving the correlation 
are given. 

2. AIRCRAFT DESCRIPTION 


The flight tests were performed on an Aerospatiale Gazelle 
SA349/2 helicopter (Fig. 1) dedicated to experimental research. The 
SA349/2 belongs to the Gazelle family of helicopters, but differs some- 
what from the 342 line of aircraft. The SA349/2 includes a comprehen- 
sive on-board data acquisition system. Table 1 provides general char- 
acteristics of the helicopter. For this flight test program, the 
aircraft was fitted with the Aerospatiale advanced geometry rotor. The 
rotor includes the classical articulated Non Articule en Trainee (NAT) 
hub (Fig. 2) and three research Grande Vitesse (GV) blades. The NAT 
hub, which is the standard Gazelle hub, is fully articulated and has a 
lag damper with nonlinear characteristics (stiffness and damping). 

The GV blades were developed to achieve a better understanding of 
rotor aerodynamics and dynamics. Features of the blades are (Fig. 3) 
0A209 advanced airfoils (9% thickness), a removable tip, and an adjust- 
able blade section c.g. location. For this flight test, a rectangular 
tip and a blade c.g. location of 24. 1 % chord were used. The 0A209 
airfoil is discussed in Ref. 2. The twist distribution (Fig. 3b) is 
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constant at the tip, linear on the main part of the blade, and varies 
11.5 degrees from the root to the first airfoiled section. 

The tail rotor for the SA349/2 is a Fenestron. The Fenestron is a 
shrouded 13-bladed rotor with one control input (collective pitch 
control) . 

Detailed technical data on the aircraft and rotors are available 
from Ref. 3. 

3. FLIGHT TEST INSTRUMENTATION 


The rotating frame instrumentation included both blade pressure 
transducers and strain gages. Strain gages for measuring pitch-link 
loads and flap and lag angle were also available. 

The blade pressure data analyzed in this paper were obtained from 
51 pressure transducers distributed between two blades, as shown in 
Fig. 4. Three groups of 20 chordwise transducers (upper and lower 
surface) were located at 75?, 88$, and 97? radial stations. The number 
of functional transducers were 18, 15, and 18 for 75%, 88?, and 97? 
radial station, respectively. 

The blade structural data in this paper were obtained from 15 
strain gages, also shown in Fig. 4. The gages (four flap-bending, seven 
lag-bending, and four torsion gages) are located at eight radial sta- 
tions on two blades. Data from additional blade strain gage instrumenta- 
tion installed on these two blades were not recorded because of slip- 
ring limitations. 

Instrumentation in the fixed frame included strain gages on the 
primary rotor control system servos and main gear box struts. Pilot and 
copilot seat vertical acceleration were also measured. The following 
flight condition parameters were recorded: ground pressure and tempera- 

ture, altitude, velocity, static temperature (in flight), static pres- 
sure (in flight), collective pitch stick position, lateral and longi- 
tudinal cyclic pitch stick position, load factor, rotor speed, sideslip 
angle, aircraft pitch, roll and yaw speed, engine power and fuel 
consumption. 

4. DATA REDUCTION 


A total of 18 consecutive rotor revolutions were required to 
obtain all of the blade pressure data. Data from each of the three 
groups of 20 pressure transducers were acquired as time-histories over a 
period of six consecutive rotor revolutions (starting from 0 degrees 
azimuth). The pressure signals were digitized at a frequency of 
180 samples/rotor revolution. The data were then converted into physi- 
cal units and transformed into aerodynamic coefficients. The strain 
gage signals were obtained as harmonic analysis coefficients (10 harmon- 
ics) calculated for one rotor revolution. The blade strain gage data 
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were acquired during one of the 18 revolutions of pressure data acquisi- 
tion. The blade strain gage data were next transformed into blade- 
bending and torsional moments after accounting for the blade weight 
contribution and flap/lag/torsion interaction components. Figure 5 
illustrates the data reduction procedure. 

5. FLIGHT CONDITIONS 


Various operating conditions were included in the orignal flight 
test envelope: (1) hover out-of-ground-effect; (2) level flights at 

305 m, 1524 m, and 3048 m (1000 ft, 5000 ft, and 10000 ft) with speeds 
from 100 km/hr to 300 km/hr; (3) steady turns at 305 m, 1524 m, and 
3048 m (1000 ft, 5000 ft, and 10000 ft) corresponding to a range of load 
factor from 1.4 to 2.0; (4) autorotation at 200 km/hr; and (5) high 
speed dive at 350 km/hr. 

From this flight envelope, 16 flight points were chosen for docu- 
mentation purposes. The selected points provide a large range of 
advance ratios (0.13 < y < 0.36) and thrust-to-solidity ratios 
(0.062 < C T /a < 0.133), as shown in Fig. 6. The data are documented in 
Ref. 3; in addition, further details of the flight test instrumentation 
and data reduction process are explained in Ref. 3. 

6. DESCRIPTION OF THE COMPREHENSIVE R0T0RCRAFT MODEL 


CAMRAD is an analysis designed to calculate rotor performance, 
aerodynamic and structural loads, aircraft vibration and gust response, 
flight dynamics and handling qualities, and aeroelastic stability. The 
analysis development is discussed in Refs. 4-6. Summaries of the vari- 
ous structural and aerodynamic models used in CAMRAD are described in 
Ref. 7. A general two-rotor aircraft with either articulated, hinge- 
less, gimballed, or teetering rotors can be modeled using CAMRAD. The 
rotor structural model is based on engineering beam theory for rotating 
wings with large pitch and pretwist. The rotor aerodynamic model is 
based on lifting-line theory, and uses two-dimensional airfoil charac- 
teristics and a vortex wake. Three levels of aerodynamic analysis exist 
in CAMRAD: uniform inflow (linear variation of inflow over the rotor 

disk), nonuniform inflow with prescribed-wake geometry, and nonuniform 
inflow with a free-wake geometry. Figure 7 outlines the aircraft trim 
solution procedure in CAMRAD. 

7. APPROACH 


The strategy for the correlation effort was to first develop a 
baseline model of the SA349/2 aircraft and rotor for CAMRAD. Much time 
and effort was spent obtaining detailed information about the GV blades, 
NAT hub, and aircraft geometry. Also, definitions, sign conventions, 
and reference systems for the data and theory were constantly checked 
for consistency. 
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Preliminary comparisons of theory with data were then per- 
formed. CAMRAD predictions were performed for seven flight conditions; 
the flight parameters for these conditions are shown in Table 2. 

Detailed results for the first three of the seven conditions will be 
presented here. These three flight points represent forward flight 
extremes: low speed, low thrust; high speed, low thrust; and high 

speed, high thrust. As a result of the preliminary CAMRAD calculations, 
several approximations and assumptions about both theory and data were 
necessary. The results presented in this paper are qualified by the 
following : 

1. Measured mean values of blade loads proved to be unrealistic; 
therefore, all strain gage time-history data presented in. this paper 
have the mean value removed. 

2. Measured blade flap and lag angles are not presented since 
their accuracy is suspect, especially at high forward speed. 

3. Measured lift coefficient is actually the normal force coef- 
ficient C n (the force perpendicular to the chord line). At small blade 
angle of attack, the difference between C a and C n is insignificant. 

For flight conditions where the blade experiences large angles of inci- 
dence, the error that is caused by the difference between and C n 
will be a small part of the overall error in predicting C^. Large 
angles of incidence occur when the blade experiences dynamic stall, and 
as will be shown, this phenomenon is not well-predicted. 

4. Measured blade moment coefficient values are unreliable, 
probably because of an insufficient number of chordwise pressure mea- 
surements . 


5. The steadiness of Condition 3 data is questionable. Varia- 
tions in the flight parameters were noticed during a period of 18 rotor 
revolutions of static data acquisition. 

6. Fuselage aerodynamic characteristics used in CAMRAD were 
obtained from a wind tunnel test of a one-seventh-scale model of the 
SA349/2 fuselage. 

7. The same aircraft c.g. location was used for all flight 
conditions. 

8. CAMRAD normally calculates the blade-bending loads at a given 
radial station, r, by computing the difference between the total aerody- 
namic and inertial moments on the blade surface outboard of r. Using 
this method for the advanced-geometry rotor, the flapwise bending 
moments were found to be very sensitive to the blade mass distribu- 
tion. This sensitivity is numerical, not physical. The blade-bending 
moments were, therefore, calculated directly from the blade curvature 
and stiffness. This method significantly improved the correlation of 
the blade loads; however, this method is not accurate where discontinui- 
ties in the blade stiffness distribution occur because of the modal 
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structural model in CAMRAD. For the GV blades, jumps in the flapwise 
and edgewise stiffness occur near the 1 0 % and radial stations. Poor 
correlation for bending moments near these radial stations is therefore 
expected . 

9. The measured control inputs were used as initial control 
settings in the free-flight trim analysis in CAMRAD. In the trim itera- 
tion, the control settings were adjusted until force and moment equilib- 
rium on the aircraft were achieved. For Condition 3, however, the 
aircraft was also trimmed to the measured power (descent). When trim- 
ming Condition 3 to force and moment equilibrium alone (level flight), 
the trimmed power level was unrealistically high. 

Baseline modeling conditions were next developed for the seven 
flight conditions. Table 3 provides baseline values used in CAMRAD. 
Modeling parameters in CAMRAD were then varied to assess the sensitivity 
on the correlation; these parameters included: number of blade-bending 

and torsion modes; wake geometry; tip-vortex core size; lifting-surface 
correction; and blade stall model. 

8. RESULTS 


Prediction of Control Inputs and Power Required 

Figure 8a presents predicted collective pitch as a function of 
advance ratio. Except for the turning flight conditions (Conditions 3 
and 7), the collective is overpredicted by approximately 0.5 degrees for 
all conditions. The nearly constant offset in collective results from a 
tens ion- tors ion coupling caused by the jump in blade twist at the 
inboard radial station (Fig. 3). This steep jump in twist was discov- 
ered to be the source of a steady, nose-down blade pitch. CAMRAD calcu- 
lates a pitch increment from the tension-torsion coupling associated 
with the pretwist. A small positive pitch increment results from the 
negative twist inherent over most of the GV blade, and a large negative 
pitch increment results from the large positive pitch change at the 
root. The net elastic pitch change is negative and occurs at the blade 
root. The pitch increment is also steady because the change is produced 
by the centrifugal force. Since the analysis trims the helicopter to 
equilibrium flight, the primary effect is on the main rotor collective 
pitch. 


CAMRAD models this tension-torsion coupling similar to Ref. 8. 
Reference 9 presents a more general theory to account for this effect. 
Because of lack of mean torsion moment and pitch link load data, the 
predicted elastic pitch increment could not be checked directly. Fig- 
ure 8a indicates that the magnitude of the elastic pitch increment is 
slightly overestimated. Calculations incorporating the general theory 
of Ref. 9 were performed by introducing an effective blade section 
radius of gyration parameter into CAMRAD. These calculations are also 
shown in Fig. 8a. The overall effect of implementing the general theory 
is to reduce the influence of the blade pretwist in the calculation of 
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the collective pitch. The magnitude of the influence depends on the 
blade shear modulus and the useful load-bearing, cross-sectional area of 
the blade, which are used to compute the effective radius of gyration 
parameter. Because the GV blade is a composite structure, the approxi- 
mation used for the shear modulus was very rough. 

Figure 8b and 8c present lateral and longitudinal cyclic pitch, 
respectively, as a function of advance ratio. The lateral cyclic pitch 
is overpredicted by an average of less than 0.5 degrees. The calculated 
longitudinal pitch is predicted to within 2 degrees (excluding Condi- 
tion 7) . 

The prediction of total power required is shown in Fig. 9a. 
Condition 3 is not shown since measured power was used to trim the 
aircraft. Except for Condition 7, predictions for the remaining condi- 
tions are within an average of 5 % of measured values. Aircraft attitude 
is not well-predicted for most of the flight conditions (Fig. 9b), 
probably because of an inaccurate representation of the fuselage 
aerodynamics. 

8.2 Condition 1: low speed, low thrust (y a 0.14, Cj/a = 0.065) 

At low forward speed, the rotor wake tends to remain in the 
vicinity of the blades and strongly influences the blade aerodynamic 
environment. A detailed wake geometry is needed therefore for good 
aerodynamic load correlation. 

Figure 10 presents C^ as a function of azimuth for r/R = 0.75, 
0.88, and 0.97. Uniform inflow is not able to predict the details of 
the lift variation. The prescribed wake results are better, but still 
do not predict the details of the distribution. The free-wake geometry, 
however, is able to predict the trend and magnitude of the lift. At 
this low speed, the rotor wake is slowly convected away from the rotor 
allowing time for interaction between the blade and shed tip-vortices. 
This interaction is evident in Fig. 10 for r/R □ 0.97 on the 
retreating side of the rotor disk (note sharp peaks in C^ near 
270 degrees azimuth). Figure 11 shows the effect of varying tip-vortex 
core size on C^. A ratio of core size to blade chord of 0.525 was 
found to be optimum. The core size used in CAMRAD includes other 
effects of blade/vortex interaction which are not modeled by the theory; 
therefore, the core size in CAMRAD is larger than the physical tip- 
vortex core size. 

For cases where the blade passes close to a tip-vortex (as in 
Condition 1), lifting-line theory usually cannot predict blade loads 
accurately, especially at the tip. A lifting-surface solution is 
required for these cases. CAMRAD uses a correction based on a linear, 
lifting-surface theory to adjust the solution for the vortex-induced 
loading on the blade which is calculated by lifting-line theory. 

Because of computation time constraints, this correction is only applied 
when the blade and vortex are a specified distance apart. Figure 12 


7 


shows C a predictions using the lifting-surface theory correction. 

When using this correction, a more realistic tip-vortex core size is 
allowed (r c /c = 0.10 was used for this figure). For this flight condi- 
tion, the calculated correction to lifting-line theory is small. The 
effect of tip-vortex core size is greater than the effect of the 
lifting-surface correction on C^. 

The flapwise bending moment is well-predicted at the 54$ radial 
station, but not at the 20$ and 80$ stations (Fig. 13). Although the 
free-wake geometry results are a significant improvement over the uni- 
form inflow results, the free wake provides only a small improvement in 
the predictions as compared to the prescribed wake results. The edge- 
wise bending moment is predicted both in magnitude and trend for 
r/R = 0.20, predicted in trend for r/R □ 0.54, and not well-predicted 
for r/R = 0.80 (Fig. 14). The difference between wake geometries 
becomes more significant near the blade outboard regions. The torsional 
moments are not well-predicted for any radial station, as shown in 
Fig. 15. If the one-half peak-to-peak values for data and theory were 
compared, however, the correlation would be quite reasonable. The data 
show a 5/rev - 6/rev influence that is not accounted for by the 
theory. CAMRAD calculates the first torsion frequency to be 3. 9/rev. 

The predicted pitch-link loads show a large 1/rev influence for 
both uniform and nonuniform inflow, whereas the data do not (Fig. 16). 

A possible explanation is the value used for the control system stiff- 
ness in CAMRAD (a very large value was used). A better estimation of 
the stiffness may improve the correlation. 

8.3 Condition 2: high speed, low thrust (y = 0.36, C T /a = 0.071) 

At high forward speed, the rotor wake is rapidly convected down- 
stream, and so wake effects on blade loading are smaller. Figure 17 
shows that the nonuniform inflow predictions of lift are slightly better 
than are the uniform inflow predictions. However, the difference 
between the free-wake and the prescribed-wake geometry results are 
insignificant. Henceforth, predictions for this flight condition are 
based on nonuniform inflow with a prescribed wake geometry. 

The flapwise bending moments are reasonably well-predicted; 
better correlation is achieved at the outboard radial stations than at 
the inboard radial stations (Fig. 18). Using six blade-bending modes 
provides better predictions at the outboard radial stations than using 
five blade modes. The difference between using five and six blade modes 
is noticeable; however, there is no difference in results from using 
six, seven, or eight bending modes. The seventh and eighth bending 
modes are dominated by edgewise motion; therefore the flapwise moment 
predictions are not affected by these modes. The edgewise moments are 
not well-predicted on the advancing side of the rotor disk (Fig. 19), 
probably reflecting compressibility effects for this high-speed case. 

For the retreating side of the disk, the moments are better predicted at 
the inboard radial stations. Using five or six blade modes appears to 
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provide the closest correlation. Figure 20 shows that pitch-link loads 
are badly predicted. Varying the number of blade modes does not improve 
the correlation. 

Reasonable correlation is obtained for torsional moments at 20% 
and 54% radial stations (Fig. 21). Increasing the number of torsion 
modes from two to four does not improve the predictions. The 5/rev - 
6/rev influence present in Condition 1 is absent here. 

8.4 CAMRAD/Finite-difference predictions of surface pressures 

With an advancing tip Mach number of 0.87, chordwise surface 
pressure data for this flight condition reveal shocks forming near the 
blade tip. Reference 10 describes a method for predicting transonic 
flows on rotors in forward flight. The predictions are made via a loose 
coupling (file transfer) between CAMRAD and a small-disturbance, finite- 
difference code (FDR). Recently, comparisons were made of results from 
FDR and several other finite-difference codes coupled with CAMRAD 
(Ref. 1). Results obtained by using one of these codes, a full- 
potential rotor code (FPR), in addition to results from FDR, will be 
presented here. FPR is described fully in Ref. 11. The code solves a 
three-dimensional, conservative, full-potential equation to calculate 
unsteady rotor flows. 

As explained in Ref. 10, CAMRAD first computes the blade partial 
angle of attack for the advancing side of the disk. The finite- 
difference code then computes the blade lift (for r/R > 0.5) for the 
advancing side of the disk. This is defined here as the end of the 
first iteration. The C^'s computed by the finite-difference code are 
then used by CAMRAD to compute new blade angles, which are used by the 
finite-difference code to compute a new C^ distribution (end of second 
iteration). The iteration continues until the difference between the 
C^ calculated by CAMRAD and the finite-difference code becomes small. 

Results from Ref. 1 are for rigid blade motion. As shown in 
Table 3, six bending modes and two torsion modes were used for the 
calculations presented here in Figs. 22 to 25. 

Figures 22 and 23 present lift coefficient predictions for the 
advancing side of the rotor disk for FPR and FDR, respectively, for 
Condition 2. As stated earlier, blade lift information is passed from 
the finite-difference code to CAMRAD for the advancing side of the disk 
only. For the retreating side of the disk, lift information from the 
previous CAMRAD iteration is used. The step change in lift at 0 and 180 
degrees azimuth can excite unwanted blade motion, especially for an 
elastic blade. In Fig. 22 and 23, the second iteration predictions are 
slightly oscillatory in nature. To avoid this problem, the step change 
in lift at 0 and 180 degrees azimuth should be made small. Ideally, 
inflow and lift information for the entire rotor disk should be 
exchanged between CAMRAD and the finite-difference code. 
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Figures 22 and 23 show that using the coupling procedure between 
CAMRAD and either finite-difference code does not improve the prediction 
of Cl (for this flight condition). In fact, the lift is underpredicted 
at all locations, except in the region of = 120 to 160 degrees 
azimuth for FPR and in the region ii> = 100 to 140 degrees azimuth for 
FDR. In these two regions, the second iteration results overpredict the 
lift. Figure 23 shows that FDR significantly underpredicts the lift 
near 0 degrees and 1 80 degrees azimuth for all three radial stations. 

In general, lifting-line theory does an adequate job of predicting the 
section lift. The reason for underpredicting the lift in Figs. 22-23 
with both finite-difference codes is difficult to determine without a 
tip-path-plane angle-of-attack measurement. 

Figure 24 and 25 show chordwise surface pressure coefficient (Cp) 
predictions for the three radial stations at three azimuth locations for 
FPR and FDR, respectively. The pressure predictions are again for 
Condition 2. As shown in Fig. 24, the shape of the Cp distribution 
predicted by FPR compares well with the data for all locations. The 
upper surface pressure peak is consistently underpredicted (because the 
calculated lift is too low). There is little difference between the 
first and second iteration results except at i|> = 150 degrees 
(reflecting the trend in lift shown in Fig. 22). At i|> = 150 degrees, 
the upper surface pressure peak is closely matched using the second 
iteration results (although the lower peak is better predicted using the 
first iteration results). Figure 25 shows similar results for the FDR 
predictions. The shape and magnitude of the pressure distributions are 
not as well predicted as compared to the FPR results. In particular, 

FPR predicts the lower surface pressure peak better than FDR. As in 
Fig. 24, differences between the first and second iteration results are 
small. FPR reasonably predicts the shock locations for r/R = 0.88 and 
0.97 at ♦ : 150 degrees (Fig. 24 b) and c)). FDR, however, does not 
predict these shock locations (Fig. 25 b) and c)). 

8.5 Condition 3: high speed, high thrust (y = 0.34, C-j./o = 0.125) 

Although not shown here, there is little difference between the 
prescribed- and free-wake results for this condition. Since the 
prescribed-wake analysis is computationally less expensive than the 
free-wake analysis, calculations were performed using the prescribed- 
wake geometry. The high thrust level of this flight condition leads to 
stall on the retreating side of the rotor disk (especially near the 
blade inboard region). The rapid blade pitch motion causes a delay in 
the occurrence of separated flow. Once the blade has reached stall, a 
vortex is shed from the leading edge resulting in large transient blade 
loads. Therefore, a dynamic stall model is used in the calculations for 
this flight condition. Several models of dynamic stall are available in 
CAMRAD. In the model used here (based on Refs. 12-14), a correction 
that is proportional to a is applied to a before the blade section 
lift, the drag, and the moment coefficients are evaluated. This model 
accounts for the hysteresis of dynamic stall, but not for the leading- 
edge vortex effects. 
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For the T5% radial station, Fig. 26 shows that the static stall 
model is not able to predict the complex behavior on the retreating 

side of the rotor disk. The dynamic stall model provides a better 
prediction in this region of the disk. The difference between the two 
stall model predictions on the retreating side of the disk becomes less 
significant for the outer radial stations since the blade becomes less 
stalled near the blade tip. As an additional comparison, one revolution 
of unaveraged data is included. Attempts were made to improve the 
predictions by varying the correction term to the blade angle of attack; 
however, significant improvements were not obtained. A more sophisti- 
cated stall model will be necessary to predict the detailed behavior 
of C a . 


Figure 27 shows that the flapwise moments are reasonably well- 
predicted. No significant improvement is gained by using the dynamic 
stall model instead of the static stall model. On the other hand, the 
edgewise moment at 20% radial station is better predicted by the dynamic 
stall model on the retreating side of the disk (Fig. 28a), where the 
blade is most stalled. The edgewise moments at the outer radial sta- 
tions are not well-predicted by either stall model. For all three 
radial stations shown in Fig. 29, the torsional moments are not well- 
predicted on the retreating side of the disk by either model. Figure 30 
again shows the poor correlation of pitch-link loads. Using the dynamic 
stall model does not improve the correlation. 

9. RECOMMENDATIONS FOR CORRELATION IMPROVEMENT 


1. Since the one-seventh-scale data does not represent the full- 
scale fuselage aerodynamic characteristics exactly, additional calcula- 
tions were performed to assess the sensitivity of the correlation to 
fuselage drag. For Condition 2, changing the fuselage drag had a pro- 
nounced effect on the trimmed control inputs, power required, and blade 
loads. The predictions were also somewhat sensitive to the fuselage 
moment characteristics. For future correlations with this data set, 
using full-scale fuselage aerodynamic characteristics would improve the 
accuracy of the predictions. 

2. Since the blade flap angle measurements were suspect, the 
rotor position relative to the fuselage was unknown. If CAMRAD were 
trimmed to the measured tip-path-plane angle, the resulting predictions 
could improve. 

3. Incorporating a more sophisticated dynamic stall model into 
CAMRAD may improve the predictions of C^ for stalled regions of the 
blade. 

10. CONCLUSIONS 


A summary of the correlation results for aircraft performance and 
trim and the three flight conditions are listed below. 


10.1 Aircraft performance and Trim 


1. CAMRAD is able to predict aircraft performance to within an 
average of 5 % for the level flight (low thrust) cases; performance for 
the turning flight (high thrust) cases is not well-predicted. 

2. The collective pitch is consistently overpredicted by about 

0.5 degrees. This overprediction is caused by the steep jump in twist 
at the blade inboard section which produces a steady, nose-down blade 
pitch increment. Introducing an effective blade section radius of 
gyration parameter reduces the effect of the pretwist on the calculated 
collective . 

3. Lateral cyclic pitch is overpredicted by roughly 0.5 degrees 
and longitudinal cyclic pitch is predicted to within 2 degrees (both 
excluding Condition 7). Aircraft attitude is not well-predicted. 

10.2 Condition 1 (y = 0.14, Cp/o = 0.065) 

1 . Blade tip-vortex interactions dominate the wake effects on 
blade loads. Tip-vortex core size has pronounced effects on the C^ 
predictions. 

2. The free-wake geometry is able to model the blade tip-vortex 
interactions. 

3. Use of a lifting-surface theory correction has little effect 
on the C^ prediction. 

4. Nonuniform inflow is required to predict blade structural 
loads. The free wake is not important in predicting flapwise moments, 
but is required for calculating edgewise bending moments at outer radial 
stations. Torsional moments and pitch-link loads are not well- 
predicted . 

10.3 Condition 2 (y = 0.36, Cp/o = 0.071) 

1. Nonuniform inflow is required for reasonable predictions of 
blade aerodynamic and structural loads; however, there is little or no 
difference between the free- and prescribed-wake results. 

2. Use of five or six blade-bending modes provide the best 
predictions of bending moments. Flapwise moments were well-predicted 
for all radial stations. Edgewise moments are not well-predicted on the 
advancing side of the rotor disk; better correlation is achieved on the 
retreating side of the disk at inboard stations. 

3. Variation of the number of blade-bending modes had no effect 
on the poor prediction of pitch-link loads. 

4. Torsional moments were reasonably predicted at the inboard 
radial station. Increasing the number of torsional modes had a small 
effect on the predictions. 
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5. Coupling CAMRAD with either FPR or FDR did not provide a 
better prediction of blade lift as compared to lifting-line theory 
predictions. 

6. In general, the FPR code provided better predictions of blade 
surface pressures than the FDR code. The shape of the pressure distri- 
butions were well-predicted by FPR. However, both FPR and FDR underpre- 
dicted the magnitude of the upper surface pressures. 

7. Shock locations were reasonably well-predicted by FPR. 

10.4 Condition 3 (v = 0.34, C^./a = 0.125) 

1 . The difference between prescribed-wake and free-wake results 
is small. 

2. A dynamic stall model provides a better prediction of C^ on 
the retreating side of the rotor disk for the blade inboard sections 
than a static stall model. The detailed behavior of C^, however, is 
not well-predicted by using the dynamic stall model. 

3. The flapwise moments are reasonably well-predicted, although 
the dynamic stall model does not offer significant improvements over the 
static model results. The edgewise moments at the inboard radial sta- 
tions are more sensitive to the stall model used; the edgewise moments 
at the outboard radial stations are not well-predicted. 

4. Torsional moments on the retreating side of the disk and 
pitch-link loads are not well-predicted by either stall model. 

11. VALUE OF PERFORMING CORRELATION STUDIES 

The benefits of performing extensive correlations, such as pre- 
sented here, can be assessed in several ways. 

1. The theoretical model is checked for validity over a wide 
range of flight conditions. 

2. Limitations of the analysis are identified. 

3. The data set is validated somewhat by the analysis. Data 
reliability can be verified by rechecking instrumentation, data acquisi- 
tion systems, and data reduction processes; however, gross errors in the 
data are easily identified by performing correlations with analyses. 
Errors made consistently throughout the data reduction process can also 
be identified by thorough correlation with a validated analysis. 
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TABLE 1. SA349/2 AIRCRAFT CHARACTERISTICS 


Aircraft gross weight 

2000 kg 

Main rotor 

number of blades 

3 

radius 

5.25 m 

blade chord 

0.35 m 

solidity 

0.0637 

rotational speed 

387 rpm 

Tail rotor rotational speed 

5919 rpm 

Engine type 

Astazou XIV H rated at 550 kw 
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TABLE 3. BASELINE CAMRAD CONDITIONS 


Condition 

Tip vortex 
core size, 

r c /c 

Number of 
bending 
modes 

Number of 
torsion 
modes 

Stall model, 
MR/TR 

Wake 

geometry 

1 

0.525 

6 

2 

static/static 

free 

2 

0.525 

6 

2 

static/static 

prescribed 

3 a 

0.525 

6 

2 

dynamic/none 

prescribed 

4 

0.525 

6 

2 

static/static 

prescribed 

5 

0.525 

6 

2 

static/static 

prescribed 

6 

0.525 

6 

2 

static/static 

prescribed 

7 

0.525 

6 

2 

dynamic/none 

prescribed 


a Trimmed to measured power level. 



Fig. 1. SA349/2 Helicopter. 


r 






Fig. 3. Grande Vitesse (GV) blade characteristics, (a) 0A209 airfoil 
profile, (b) Blade twist distribution. 
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SA 349 2 
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Fig. 5. Date reduction process. 












ADVANCE RATIO 

Fig. 6. Thrust and speed range of flight conditions. 
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LEVEL FLIGHTS 1 TURNING FLIGHTS 


TRIM 


UNIFORM INFLOW 



(a) Inflow Analysis Levels 

TRIMMED SOLUTION 

ITERATE CONTROLS TO TRIM 
PERIODIC MOTION AND AIRLOADS 


CIRCULATION ITERATION 
(INDUCED VELOCITY FROM LOADING) 

MOTION ITERATION 
ROTOR MOTION 

AZIMUTHAL STEP 

RADIAL INTEGRATION 
SECTION AERODYNAMICS 
UPDATE MOTION HARMONICS 
TOTAL ROTOR LOADS 


BODY VIBRATION 
TEST MOTION CONVERGENCE 
TEST CIRCULATION CONVERGENCE 


TEST TRIM CONVERGENCE 


(b) Trim, Circulation, and Motion Iterations 

Fig, 7. CAMRAD solution procedure. 
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□ MEASURED 
• CALCULATED 

A CALCULATED (BASED ON REF. 9) 


O 

O 

8.8 





AIRCRAFT PITCH, deg POWER, hp 



ADVANCE RATIO 


Fig. 9. Trimmed aircraft characteristics (numbers correspond to flight 
condition), (a) Performance, (b) Aircraft attitude. 
















UNIFORM INFLOW 



80 160 ’ 240 ’ 320 

AZIMUTH, deg 


Effect of wake geometry on pitch-link loads - Condition 1 
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cient predictions - Condition 2. 
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(c) r/R = 0.97: azimuth = 30, 90, 150 degrees 

Fig. 25. Concluded. 
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